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We use spatially resolved scanning tunneling spectroscopy in Na(Fei-xCox)As to investigate the 
impurity effect induced by Co dopants. The Co impurities are successfully identified, and the spatial 
distributions of local density of state at different energies around these impurities are investigated. 
It is found that the spectrum shows negligible spatial variation at different positions near the Co 
impurity, although there is a continuum of the in-gap states which lifts the zero-bias conductance 
to a finite value. Our results put constraints on the S=b and S++ models and sharpen the debate 
on the role of scattering potentials induced by the Co dopants. 
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I. INTRODUCTION 

According to the Anderson's theorem 1,2 , the Cooper 
pairs with singlet pairing can survive in the presence of 
non-magnetic impurities, while the magnetic impurities 
are detrimental to superconductivity yielding a unique 
pattern of the electronic local density of states (LDOS). 
This view has been well established in the conventional 
superconductors 3 . In superconductors with gap nodes, 
even the nonmagnetic impurities can be strong pair 
breakers 2,4 . In the iron-based superconductors 5 , the 
pairing mechanism remains as a hot and unresolved is- 
sue. It is proposed that the pairing is mediated by 
the antiferromagnetic (AF) spin fluctuations resulting 
in a S± pairing 6,7 , which now becomes a widely per- 
ceived picture. The experimental evidence that di- 
rectly supports S± pairing is still quite rare. One 
was drawn from the scanning tunnelling spectroscopy 
(STS) and the quasi-particle interference measurements 
in Fe(Se,Te) 8 . Another indirect evidence is the obser- 
vation of the resonance peak of the imaginary part of 
the spin susceptibility at (7r, it) in the inelastic neutron 
scattering experiment^. Theoretically it was argued that 
this type of pairing should also be fragile to the non- 
magnetic impurities^"— . One of the puzzles concern- 
ing this model is why superconductivity survives up to 
a relatively high temperature in some doped samples, 
such as the Ba(Fei_ a ,T a ,) 2 As 2 (T=Co, Ni, Ru, etc.), 
where the dopants are supposed to be the pair break- 
ers. The spatial and energy dependent LDOS induced 
by the impurity scattering can be regarded as the fin- 
gerprints for a formal check on this peculiar pairing 
modeliSr— which can be measured by the STS. After 
cleaving an iron-based superconductor, one usually en- 
counters a polar and messy surfac e - 16 ! 17 , while an atom- 



ically resolved image was sometimes observed in a re- 
constructed surface^ - — . However, it is very difficult to 
locate a well-defined quantum impurity and investigate 
the related LDOS on such surfaces with an exception 
in Ca(Fei_ a: Co a: )2As2^. In this paper, we report the 
success in locating the Co impurities and measuring the 
spatially resolved STS at different energies and tempera- 
tures obtained from the Na(Fei_ a: Co a: )As single crystals. 
Our results indicate an unexpected weak variation of the 
STS and LDOS when going through a Co impurity site. 
Combined with the theoretical calculations, we propose 
explanations for this dichotomy (finite LDOS but weak 
variation). Our findings put constraints on the theoreti- 
cal models of the pairing state. 



II. EXPERIMENTS 

High quality Na(Fei_ a: Co a: )As single crystals were syn- 
thesized by the flux method 21 . The samples were cleaved 
at room temperature in an ultra- high vacuum with a base 
pressure better than 2 x 10 -10 Torr, then transferred into 
the scanning tunneling microscopy (STM) head imme- 
diately and cooled to a desired temperature. The STS 
were measured with an ultra-high vacuum, low tempera- 
ture and high magnetic field scanning probe microscope 
USM-1300 (Unisoku Co., Ltd.). In all STM/STS mea- 
surements, Pt/Ir tips were used. The surface topogra- 
phies were recorded using a bias voltage of Vbias = 40 mV 
and tunneling current of I t = 100 pA. To reduce the 
noise of the differential conductance spectra, a lock-in 
technique with an ac modulation of 0.1 mV at 987.5 Hz 
was typically used. 



III. RESULTS 



A. Sample characterization 

Figure QJa) shows the temperature dependence of the 
volume dc magnetization (M) after zero-field-cooling 
(ZFC) and field cooling (FC) processes. The magnetiza- 
tion on the samples of two doping levels showed sharp su- 
perconducting transitions with T c « 16.8 K for x = 0.05 
and T c « 21.2 K for x = 0.025. Since the single crystals 
are thin and plate like, the demagnetizing factor is rel- 
atively large; that is the reason the ZFC magnetization 
values of these two samples are different. The error in 
the measurement of dimensions, especially the thickness, 
could give some error in the calculation of the ZFC sus- 
ceptibility. Figure 0Jb) shows the resistive transitions of 
the two samples. It is clear that the sample with more 
charge doping level (x=0.05) is lower than that at the 
optimal doping (x=0.025) although their T c values are 
opposite. The different values of the resistivity at 300 K 
may be related to the charge doping effect in the samples. 



B. Identification of Co dopant atoms in 
Na(Fei_ x Co x )As 

Figures [2ja) and EJb) present the topographic image 
of the cleaved surface of samples with (a) x = 0.025 and 
(b) x = 0.05. The background shows a well-ordered 
square lattice with a constant of 3.80 A. Because the 
Na atoms are distributed equally into two layers between 
two neighboring As-layers, the cleaving occurs between 
the two Na layers leading to an un-polar (charge neutral) 
surface 22 . In both samples, beside the observed square 
lattice, one can see some rectangular blocks with two 
unit cells constructed by six Na atoms (hereafter named 
a 2 x 1 block) which are a bit brighter than the back- 
ground as shown in Figs. [2fa) and 2(b). These blocks 
align along either [100] or [010] directions. The bright- 
ness is enhanced when many blocks overlap each other 
[marked by an arrow in Fig. Efa)]. After considering the 
arrangement of Na atoms in one 2x1 block, we argue 
that each block is corresponding to one individual Co im- 
purity atom. Two reasons can be given to support this 
argument, as addressed below. 

Firstly, when one Co atom is doped to the Fe site, 
it naturally neighbors with six Na atoms on the sur- 
face layer, forming a 2-unit-cell rectangular (2x1) block 
as illustrated by Fig. Efe), the orientation of the 2x1 
block is determined by the selective positioning of the 
Co atom on the a-axis or 6-axis of the Fe-As square 
lattice. Secondly, counting the ratio between the num- 
ber of the typical 2x1 blocks and the total number 
of Fe and Co atoms, we find that it is very close to 
the doped concentration of Co, i.e., the nominal value 
of Co/(Fe±Co) in the samples. We counted the num- 
ber of 2 x 1 blocks as shown in Fig. [2ff) and (g). The 
very bright regions on the over-doped Na(Feo.95Coo.o5)As 



E 
jj 

E 

CD 



x=0.025ZFC ; 

— — x=0.025FC 
x=0.05ZFC 

— ° — x=0.05FC 



L^ 



r Na(Fe 1 x Co x )As. 



10 15 20 25 

7(K) 




J i l i l i l_ 



50 100 150 200 250 300 
7(K) 



FIG. 1: (Color online) (a) Temperature dependence 
of volume magnetization of Na(Feo.975Coo.o2s)As and 
Na(Feo.95Coo.o5)As single crystals after zero- field cooling 
(ZFC) and field cooling (FC) at 10 Oe. The difference in 
the ZFC magnetization of the two samples comes from the 
different demagnetizing effect, (b) Temperature dependence 
of resistivity for the two Na(Fei_ cc Co ;E )As samples with dif- 
ferent doping levels. 



sample are originated from overlapping of several 2x1 
blocks which can be distinguished from the differential 
figure as shown in Fig. EJh). The stacking blocks are 
distinguished from different height of the Na atoms and 
contrast level of the experimental image. The quanti- 
ties of Co atoms counted on Na(Feo.975Coo.o25)As and 
Na(Feo.95Coo.o5)As samples are 64 and 110 respectively 
with the counting error of 10%, while the total number 
of Fe sites should be 2677 in such 14 nmxl4 nm area. 
Only one bright region with the question mark at the 
bottom edge on the Na(Feo.95Coo.o5)As sample as shown 
in Fig. [2fg) hasn't been taken into account, due to the 
uncertainty. The corresponding ratio of doped Co atoms 
over all the Fe±Co sites Co/(Fe±Co) on the sample with 
x = 0.025 « 2.39 ±0.24%; while on the sample x = 0.05, 
we get Co/(Fe±Co) w 4.11 ±0.41 %. The two calculated 
ratios are comparable with their chemical doping levels 
2.5% and 5%, respectively. This could not be achieved 
by accident. We also checked the height where several 
blocks overlap each other. As shown in Figs. [3] (a), (b), 
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FIG. 2: (Color online) (a, b) High-resolution topography in a 
region of 14 x 14 nm 2 on Na(Fei_ x Co x )As samples with x = 
0.025 (a) and x = 0.05 (b). (c, d) The spatial dependence of 
height z measured along the white lines marked in (a) and (b). 
(e) Illustration for atoms near the Co-impurity sites. It is clear 
that each Co atom neighbors with six Na atoms (highlighted 
by the yellow circles) , forming the 2x1 rectangular block (see 
text), (f, g) Illustration for counting the 2x1 blocks on (a, 
b), and the 2x1 block patterns are marked with light-blue 
rectangles, (h) The differential picture of (b) which is helpful 
to make the counting. 



(d), and (e), we present the height distributions crossing 
one single 2x1 block in two perpendicular directions. 
We used these distributions as basic cells to simulate the 
height where several blocks overlap by a simple summa- 
tion. The simulation is undertaken based on the distri- 
bution of the 2x1 blocks shown in Fig. [2ff) and (g). 
The accumulated height equals to the summation of the 
total contributions given by several blocks demonstrated 
by the frames. We test this simple method along three 
typical lines in Fig. [2ff) and (g). One can see that the 
consistency between the simulated results and the data is 
remarkable. This clearly convinces that the very bright 
spots are constructed by several 2x1 blocks overlapping 
each other, which validates the counting method we pro- 
posed here. This result may also suggest that the local 
work function on six Na atoms of the 2x1 blocks is re- 
duced, which results in a linear enhancement of the height 
when the blocks overlap each other. Such conclusion can 
be derived by the expression of tunneling current if the 
second order small quantities in the Taylor expansion are 
omitted. 

Recently, Erwin and Mazin have performed first prin- 
ciples calculations of the map of tunnelling current in 
NaFeAs with an isolated surface Co impurity 23 , and 
found that the calculated topographic image closely re- 
sembles ours and reproduces the 2x1 block that we as- 
signed to a Co impurity. These 2x1 blocks seem to 
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FIG. 3: (Color online) Simulations of the height on overlap- 
ping of 2 x 1 blocks. (a),(b),(d) and (e) are line profiles of a 
single block marked by green and violet dashed/solid lines in 
Figs.[2^f) and (g), while (c), (f) and (g) show three simulated 
line profiles of combined altitude by overlapping blocks illus- 
trating in Figs.[2{f) and (g). One can find that the simulated 
lines match well with the experimental data. 



be slightly higher (about 30 pm) than the surrounding 
atoms, as shown by the line-cut in Figs. Efa) and (b) 
with the landscape shown in Figs. [2fc) and (d), respec- 
tively. At some positions, we see dark areas characterized 
by a lower height, as shown by the line-cut in Figs. [2fb) 
and (d). These are probably induced by the vacancies of 
the atoms beneath the surface layer. 



C. Impurity state on Co-dopant site 

1. ST S far away from the Co-impurity 

In the sample with x = 0.05 we measured the evolution 
of tunnelling spectrum with temperature and showed it 
in Fig. |4^a). At low temperatures, two coherent peaks 
near the gap edges can be clearly seen. The asymmetric 
structure of the STS curve is induced by the background, 
which can be evidenced by the data measured at 25 K 
(> T c = 16.8 K) as shown in Fig. G^b). In Fig. g^c), 
we present the experimental STS at 1.7 K that normal- 
ized by the data at 25 K. On the STS curve, there is 
a finite value of zero-bias conductance, which seems to 
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FIG. 4: (Color online) (a) The STS measured on 
Na(Feo.95Coo.o5)As at a position far away from the Co- 
impurities at temperatures from 1.7 K to 55 K. (b) STS taken 
at 1.7 K and 25 K. (c) The STS measured at 1.7 K normalized 
by that of 25 K (symbols). The solid line is a fit according 
to the Dynes formula with a gap of 4.5 meV and a scattering 
rate of T = 0.7 meV. 



be a common feature of Co-doped samples 24 26 , and will 
be addressed later. No feature of two gaps is observed 
here. This is consistent with the angle-resolved photo- 
emission spectroscopy measurements on the same sam- 
ple, which indicates the same gaps, ~ 5 meV, on both 
the electron and hole pockets. 27 The normalized data 
was fitted with the Dynes model 28 with the gap func- 
tion A(0) = A[(l — a) cos 2$ + a] to fit the experimental 
data. First we try the isotropic S-wave with anisotropic 
factor a = 1, as shown in the Fig. Ef(a). The result is 
shown by the green curve. In Fig. [5(b), we show the fit- 
ting curves with different anisotropy a and T. It is found 
that an optimized fitting can be obtained when using 
a = 0.65, i.e., A(0) = 4.5(0.35 cos2<9 + 0.65) meV and 
scattering rate T = 0.7 meV. One point we should men- 
tion is that, there is a slight enhancement of dl/dV at 
about -5 meV even above T c , e.g., at 25 K. This feature 
is smeared out at elevated temperatures and vanishes at 
about 55 K. Interestingly this feature seems to be linked 
with the asymmetry of the STS, for example, the coher- 
ence peak at —A = —4.5 meV is higher than that at 
+A = 4.5 meV, even at 1.7 K. The relationship between 
the high temperature feature and the asymmetric super- 
conducting coherence peaks remains to be understood. 



2. STS across a Co-impurity 

In order to have a detailed comparison between the 
LDOS patterns as expected by the theoretical models and 
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FIG. 5: (Color online) STS measured on Na(Feo.95Coo.os)As 
at 1.7 K and the fitting with different gap functions (a) and 
different values of fitting parameters for anisotropic s-wave 
gaps (b). 



that induced by the Co-impurity, we select a 2 x 1 block 
and do the spatial resolved STS measurements with spac- 
ing of every half lattice constant, as marked by the dots 
in Fig. [6(a). The raw data of STS are shown in Figs. [6(b) 
and [6(c) for the two perpendicular directions. The data 
measured right on top of the Co-impurity is shown here 
by a red curve. It is astonishing that the shape of the 
STS does not change in any noticeable way upon moving 
through the Co-impurity. The value of dl/dV taken from 
these STS curves at three typical voltages (Vbias = —4.3, 
0, and 4.3 mV) are shown in Figs. [6(d) and [6(e) with 
respect to the positions. One can clearly see that the 
LDOS at either zero energy or the gap edges are quite 
uniform. 

In order to examine the spatial dependence of the spec- 
trum on the Co site and nearby, we show the spectra at 
different positions in Fig. [7(a) and [7(b). We also mea- 
sured the STS just on the Co site in a large bias voltage 
range as shown in Fig. [7(c). One can find an asymmetric 
background which is similar to the spectra measured on 
LiFeAs 22 . The theoretical simulation indicates that there 
are clear Co dopant induced resonances at -800 mV and 
200 mV in BaFe2As2^, and some special characteristics 
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FIG. 6: (Color online) (a) Illustration of the atomic struc- 
ture based on the topographic image (shown by the inset) of 
Na(Feo.95Coo.o5)As at 1.7 K, with a Co-impurity highlighted 
by a red circle surrounded by Na atoms (big yellow and blue 
circles), (b, c) Spatially resolved STS measured along the two 
arrowed lines in (a) with a spacing of every half lattice con- 
stant. The red lines show data measured just at the center 
of the 2x1 rectangular block, (d, e) The value of dl/dV 
determined near the gap edge (-4.3 meV and 4.3 meV) and 
zero bias when crossing the Co-impurity point. 



at about 150 mV was observed and suggested to image Co 
atoms on the surface of slightly Co-doped CaFe2As2 20 . 
However we didn't find any characteristic features within 
±200 meV at either the Co site, or somewhere far away 
from the Co sites. To investigate the spectrum affected 
by the black holes and the very bright spots originated 
from the overlapping of several 2x1 blocks in Fig. Efb), 
we measured the spectra on such places and showed them 
in Fig. E{d). One can see that the gap values as well 
as the bottom heights are almost the same at different 
positions. The superconducting coherence peaks on the 
spectrum measured on a bright spot are suppressed com- 
paring to the one measured at the point without any fea- 
tures. The black holes seem to affect the spectrum very 
little which suggests that it arises from the depletion of 
some Na atoms on the top surface. 



3. Theoretical simulated results of the LDOS on the Co site 

In order to understand the spatial dependence of the 
spectrum on the Co site and nearby in our experiments, 
we did theoretically calculations based on the S=b and 
S++ models. First we analyze the gap functions of the 
S± wave superconductivity obtained in random phase 
approximation (RPA). We use the multi-orbital Hubbard 
model proposed by Kuroki et ali, which is down-folded 
from the first-principle calculation. In this scheme, five 
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FIG. 7: (Color online) (a, b) The tunnelling spectra normal- 
ized by the normal state one measured at different positions: 
1, on site of the impurity, 2-4 corresponds to those measured 
at the 2nd, 3rd and 4th lattice constants away in two perpen- 
dicular directions, (c) STS measured on Na(Feo.95Coo.os)As 
at 1.7 K. It is found that the STS measured just on the Co 
site is very similar to the one measured far away and it has no 
special feature in the voltage range from -200 mV to 200 mV 
besides the superconducting gap at low energies. An asym- 
metric background is seen here which may be induced by the 
multi-band and the band edge effect, (d) STS measured on a 
very bright spot, a black hole and somewhere far away from 
such features. One can see that the energy gap values are 
almost the same at these specific locations. 



kept, and the hopping integrals between these orbitals for 
the LaFeAsO (1111) system is used. This will be justified 
since the Fermi surface for NaFeAs near k z = is sim- 
ilar to that of LaFeAsO, although NaFeAs has stronger 
three-dimensionality than 1111 system. The effect of the 
z-direction hopping on the impurity scattering is a fu- 
ture problem. In this paper, we do not go into details of 
the material-dependent model, but instead we compare 
the experiment and theory in a simple model which con- 
tains essence of the multi-band unconventional (S=b and 
S++) superconductivity. It is found that, when we use 
orbital representation instead of band representation, the 
gap functions obtained in RPA can be reproduced quite 
well by using the real-space short-range pairings (up to 
next-nearest-neighbor sites). The details are explained 
elsewhere. 29 Note that, in the RPA calculation around 
T c , the absolute values of the gap functions are not deter- 
mined, and only their relative sizes are obtained. There- 
fore, we adjust the largest gap value to be consistent with 
the experiment. On the other hand, for the S++ case, 
the gap functions are obtained from the analysis of the 
five-orbital attractive Hubbard model. 



The impurity is treated as a single site in the two- 
dimensional square lattice with an on-site impurity po- 
tential. The orbital dependence of the potential is ne- 
glected. The value of the impurity potential is chosen 
as I — 1.0 eV— and / = —0.35 eV 14 according to the 
previous first-principle calculations. The gap functions 
around the impurity site are assumed to be spatially uni- 
form for simplicity. The effect of this simplification is 
discussed shortly. The actual calculation of LDOS is per- 
formed using Chebyshev polynomial expansion method 30 
with 900 x 900 lattice sites and 3600 polynomials. The 
Lorentzian-type broadening kernel is used and the broad- 
ening factor is determined so as to make 7/A max consis- 
tent with experiment. Because of this broadening, LDOS 
at E = remains finite even in the S++ case, which 
comes mainly from the tail of the peak at the gap edge. 
The electron number is fixed at n = 6.05 by adjusting 
the chemical potential. 

There are two important factors that should be in- 
cluded when we compare the theoretical calculations with 
experimental data. First, there is a band renormaliza- 
tion effect due to electron-electron correlation which is 
neglected in the first-principle calculations. In iron pnic- 
tides, the effective mass obtained in ARPES is larger 
than that obtained in the first-principle calculation by 
a factor of about 2.5. Therefore, in order to take ac- 
count of the band renormalization factor (not included 
in the first-principle calculations), the hopping integrals 
and site-potentials in the model Hamiltonian should be 
reduced by a factor 2.5. By considering this effect, we 
reduce the energy scale of the obtained results by factor 
2.5, for simplicity. The other important factor is that 
the experimental tunneling spectra were taken under the 
condition of constant tunneling current. In all the ex- 
periment measurements, the spectra were taken with the 
same initial tunneling current (e.g., 100 pA) at some bias 
voltage Vo (e.g., 40 mV). So there was a normalization 
process using the integral area of sample LDOS from to 
Vq of each curve if we simply considered the tip LDOS as 
a constant. 19 Such normalization makes very little change 
of the tip height which can be omitted because the tun- 
neling current decays exponentially with increasing of 
height. This means that LDOS obtained theoretically 
should be normalized in accordance with the experimen- 
tal condition of constant current. This is carried out by 
rescaling the amplitude of LDOS for each site so that 
the energy integral of LDOS between E = to -20 meV 
becomes a constant for the negative parts of calculations 
have fewer features. 

After these treatments, LDOS shown in Figs.[8(e)-|8(g) 
are obtained. For the case of S=b and scatting potential 
I — 1.0 eV [Fig. 15(e)] and at the impurity site, the simu- 
lated data indicate a strong in-gap bound state peak at 
zero energy as shown in the previous paper—. For other 
cases, the bound state peaks are absent and shape of the 
spectrum changes when it is moving away from the im- 
purity site. From Figs. [5(f) and[8(h), we can see that the 
results with I = -0.35 eV for S++ (Fig|8(h)) are quite 
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FIG. 8: (Color online) (a-d) Un-normalized STS calculated 
in a simplified model based on S++ and S=b pairing gaps at 
several sites near the impurity with the same spatial registra- 
tion as the experiment (based on the band parameters of the 
LaFeAsO system 29 ) (e-f) Curves of (a-d) normalized by the 
integral area of LDOS from -20 mV to mV in corresponding 
to the constant-current situation in the experiment. 



consistent with the experimental data. However, we can- 
not exclude the possibilities of S=b model with / = —0.35 
eV and the S++ model with 1 = 1 eV, since the simu- 
lated results and experimental data are not very distinct. 
At the end of this part, we want to check the validity 
of neglecting the spatial dependence of the gap functions 
around the impurity site. For this purpose, we have to 
determine the gap functions around the impurity self- 
consistently by solving the Bogoliubov-de Gennes equa- 
tion iteratively just as in Ref iii. However, a large numer- 
ical cost of solving Bogoliubov-de Gennes equation limits 
the accessible energy scale. If we try to reproduce the ex- 
perimental gap size, the self-consistency calculation be- 
comes unstable or unrealistic suffered from the finite-size 
effect of the available lattice. This is just a technical rea- 
son, but we have to assume large attractive interactions 
in order to obtain self-consistent solutions. In this cal- 
culation, we solve Bogoliubov-de Gennes equation used 
in Ref. 11 that is constructed so as to reproduce the RPA 
S± results and contains attractive interactions up to the 
next-nearest-neighbor sites. On the other hand, for the 
S++ case, we use a model with short-range attractive 
interaction. We use a lattice with 28 x 28 sites with an 
impurity at the centerin order to obtain eigenvalues and 
eigenfunctions. For the LDOS calculation, we use the 
"supercell" method in which 28 x 28-site lattice with an 
impurity is treated as a "unit cell" and the whole system 
is assumed to be composed of a 13 x 13 repeat of this unit 
celliii Figures [9ja)-[9{d) show the LDOS obtained by this 
self-consistent calculation. Since the attractive interac- 
tion is large, the obtained gap function is larger than the 
experimentally observed gap values. Note that LDOS is 
not normalized and the band renormalization factor is 
not taken into account. For comparison, we show the 
LDOS without self-consistent calculation in Figs [5(a)- 
[8(d), which are identical with Figs.[8(e)-[8(h) but they are 



(a) faraway 




Energy (meV) 



Energy (meV) 



FIG. 9: (Color online) (a-d) The STS calculated in a five- 
orbital model based on S++ and S=b pairing gaps at sev- 
eral sites near the impurity with the same spatial registra- 
tion as the experiment (based on the band parameters of the 
LaFeAsO system 29 ) (e-f) Curves of (a-d) normalized by the 
integral area of LDOS from -80 mV to mV in corresponding 
to the constant-current situation in the experiment. 
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FIG. 10: (Color online) Temperature dependence of specific 
heat of Na(Feo.95Coo.o5)As. The inset shows an enlarged view 
of the same data in the low temperature limit. 



not normalized. From this comparison, we can see that 
the results do not change qualitatively even if we assume 
a site-independent gap functions [Figs [8^a)-lHJd)] com- 
pared with the full self-consistent calculation in Figsj9fa)- 
EJd). One dominant difference is a rather large peak at 
E > 50 meV in Figs^a) and E{c) (i.e., for I = 1.0 eV 
case) on the impurity site (line 1: red curve). We find 
that this large peak is due to the impurity-induced reso- 
nant state that exists even in the normal state. Although 
the impurity-induced resonant states exist in Figs. |5(a) 
and[5Jc), it is located outside of the energy-range shown 
here. 



D. Discussion on the Impurity effect by Co doping 

The experimentally measured STS seems to change 
very weakly when crossing the Co site as mentioned 
above. On one hand, we don't see the sharp bound state 
peaks at the impurity site as theoretically predicted for 
the S± pairing gap for / = 1 eV with the non-magnetic 
impurities. On the other hand, we do see a continuum 
of in-gap states, which lifts the zero-bias conductance to 
a finite value. Theoretically we find that the effect of 
Co-impurity on LDOS is similar both for S± and S++ 
case when / = —0.35 eV. One may argue that the scat- 
tering effect given by the Co- impurity is very weak i.e., 
I ~ 0. However, specific heat data shows a significant 
residual specific heat coefficient 70 = 9.8 mJ/mol-K 2 (as 
shown in Fig. [TO]) suggest that the scattering in the Co- 
doped case is significant. This phenomena occurs in other 
doped samples as well when the dopants go directly to 
the Fe sites^i. The density functional calculation of the 
local substitution of Co or Ni for Fe suggest that no extra 
charge carriers are doped into the system 32 . In this case, 
the Co-impurity should behave as an impurity scatterer 



for Cooper pairs 3 - 3 -. Provided the expectable influence 
given by the Co-impurities, it is paramount that we un- 
derstand the lack of the spatial variation of the LDOS 
spectrum nearby the Co-site within the S± model. 



Now we discuss the possible reasons for this dichotomy, 
and constraints given by our results on the theoretical 
pairing models. Although there is a discrepancy (es- 
pecially near the gap edge) between the experimental 
data and the theoretical results, the S++ model would 
give a much better consistency. However, the S=b model 
may still be relevant in the following circumstances: (1) 
The scattering potential / is around -0.35 eV, as shown 
in Fig. EJf). In the absence of precise independent as- 
sessment on the scattering potential, we assume that 
the scattering potential is relatively weak. (2) The Co- 
impurity may behave as an extended scattering center, so 
that the scattering between the electron and hole Fermi 
pockets (if would be the cause of pairing) will not be 
affected due to the finite range. Since the inter pocket 
scattering is the major cause for the pairing based on 
the S± model, that may be the reason why the super- 
conductivity survives in the presence of such impurities. 
With this picture, the obvious scattering as revealed by 
the residual resistivity and specific heat data are proba- 
bly induced by the intra-pocket scattering which can sup- 
press the superconductivity moderately, but not strongly. 
Such scattering results in the expectation of in-gap filled 
states with no sharp resonance states seen. (3) The last 
possibility would be that the pairing is induced by strong 
local super-exchange effec t 34 ! 35 which will certainly not 
be weakened so much by the non-magnetic impurity and 
relatively weak scattering. 



IV. CONCLUSIONS 

In summary, we successfully locate the doped Co 
atoms in Na(Fei_ a: Co ;c )As using scanning tunnelling mi- 
croscopy. Negligible spatial variation of the tunnelling 
spectrum near the Co impurity is observed, which is in 
sharp contrast with the theoretical simulations with the 
S± model taking 1 = 1 eV. We explain our data with the 
S++ or S± models. The S++ model can be consistent 
with the data, however, the S± model is still relevant in 
the following three cases: (1) The scattering potential I 
given by a Co dopant is relatively small, for example, I 
-0.35 eV; (2) The Co-dopants give rise to the extended 
scattering centers which will not significantly affect the 



interpocket scattering; (3) The pairing may be induced 
by local strong magnetic super-exchange effect. 
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